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Abstract 

This paper is devoted to studying obstacle avoiding patterns and cohesiveness 
of fish school. First, we introduce a model of stochastic differential equations 
(SDEs) for describing the process of fish school’s obstacle avoidance. Second, 
on the basis of the model we find obstacle avoiding patterns. Our observations 
show that there are clear four obstacle avoiding patterns, namely. Rebound, 
Pullback, Pass and Reunion, and Separation. Furthermore, the emerging pat¬ 
terns change when parameters change. Finally, we present a scientific definition 
for fish school’s cohesiveness that will be an internal property characterizing the 
strength of hsh schooling. There are then evidences that the school cohesiveness 
can be measured through obstacle avoiding patterns. 
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1. Introduction 


Fish schooling, one of animal swarming, is a commonly observed phenomenon 
that is coherently performed by integration of interactions among constituent 
fish. This remarkable phenomenon has already attracted interests of researchers 
from diverse fields includin g biology , p hysics, mathematics, computer engineer- 

li [rzi y, y, y, y , yi). 

Let us recall here some researches in the literature therein. In 2001, Ca- 
mazine et al. _0, Chapter 11 


ing (see 


results (Aoki [l[, Huth-Wissel 


raesented an idea on the basis of experimental 


14| . and Warburton-Lazarus [2^) that individual 


hsh may act following the behavioral rules: 


(a) The school has no leaders and each fish follows the same behavioral rules. 

(b) To decide where to move, each fish uses some form of weighted average of 
the position and orientation of its nearest neighbors. 

(c) There is a degree of uncertainty in the individual’s behavior that reflects 
both the imperfect information-gathering ability of a fish and the imperfect 
execution of the fish’s actions. 


Their insight is that these local rules can altogether create the coherent behavior 
of fish school. 

Vicsek et al. [27| modeled the movement as self-driven particles obeying 


some difference equations. They assumed that each individual is driven with a 
sum of an absolute velocity and an averaged velocity of nearby particles together 
with some random perturbations. Oboshi et al. [20| also modeled schooling by 
some difference equations setting a rule that each fish choose one way of action 
among four possibilities according to a distance to the closest mate. Meanwhile, 
211 and D’Orsogna et al. independently presented differential 


Olfati-Saber 


equation models, but deterministic ones, utilizing the generalized Morse function 
and attractive/repulsive potential functions, respectively. In [i^ . Gunji et al. 
considered dual interaction which produces territorial and schooling behavior. 


In 


Reynolds introduced some simple behavioral rules of animals, which are 


2 




















similar to (a), (b), (c) but deterministic ones, and are schematic rather than 
physical. 

A stochastic differential equation (SDE) model describing the process of 


schooling was presented in [26j, where we used the behavioral rules (a), (b), and 
(c) above. We then utilized the model for developing quantitative arguments 
on fish schooling in 171. 


In the real world, the environment surrounding fish school often includes 
other components such as obstacles, food resources, predators, etc. In those sit¬ 
uations, fish exhibit more complex, parallel movements such as obstacle avoid¬ 
ance, food finding, escaping from predator. It is evident that when a school 
of fish is tackled by obstacles or is hunted by predators, fish individually react 


quickly for avoiding obstacles 0 £ 
Olfati-Saber and Murray 


predators. 


22| developed the method of Reynolds [23| by 


introducing a dynamic graph of agents in the presence of multiple obstacles. 
The agents are split into several groups while approaching the obstacles. After 


passing all obstacles, they rejoin into a single group. Chang et al. [6[ introduced 
techniques of using gyroscopic forces for multi-agent systems by which the agents 
perform collision avoidance toward obstacles. A similar result has been shown 
(i.e., agents are separated into some clusters and then rejoin into a single flock). 
In the meantime, Hettiarachchi and Spears 1^ used virtual physical forces in 


composing a swarm system of robots moving toward a goal through obstacle 
fields. Robots may collide with obstacles but then they can still move toward 
the goal. 

In the meantime, a concept concerning animal swarming or grouping, namely 
cohesiveness has already been introduced since 1930s. The study during long 
years seems to show that it is not an easy problem to define a concept of the 
cohesiveness precisely and consistently. It has been conceptualized in various 
ways, but each was based on intuitive assumptions and interpretations. 


For instance, Moreno and Jennings m defined cohesiveness as the forces 


holding the individuals within the group to which they belong. French [9| noted 
that the group exists as a balance between cohesion and disruptive forces. Not 
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until 1950 was a systematic theory of group cohesiveness constructed by Fes- 
tinger et al. Q . Their definition of cohesiveness is “We shall call the total field 
of forces which act on members to remain in the group the “cohesiveness” of 
that group”. Gross and Martin Q claimed that this definition is inadequate, 
and they proposed an alternative definition as the resistance of group to dis¬ 
ruptive forces. Contemporary works almost characterize group cohesion in the 
same way (see Q). Carron B defined cohesiveness to be the adhesive prop¬ 
erty of group. Schachter et al. Q found that interpersonal attraction is the 
cement binding group members together. For the general relationship between 


0,Q,Q. 


cohesiveness and group performance, we refer the reader to 

We can however find a point of view which is common in those definitions. 
That is the bond linking group members to others and to the group as a whole. 
We believe that this common point of view may be a key feature for all inter¬ 
communicated multi-agent systems. 

The objective of the present paper is two-folds: namely, studying the fish 
schooling from a viewpoint of pattern formation of biological systems and in¬ 
troducing a scientific definition of the school cohesiveness. 

For the first objective, obstacle avoiding patterns of fish school are studied 
by newly introducing a behavioral rule for avoidance and adding its effect to our 
model in 26|. It is then observed that there are at least four obstacle avoiding 


patterns of school, i.e.. Rebound, Pullback, Pass and Reunion, and Separation 
which are performed just by tuning modeling parameters. 

For the second objective, we consider school cohesiveness as its ability to 
form and maintain the school structure against the white noises affecting the 
school. It is therefore defined as an internal nature of the school, indepen¬ 
dent of external effects. We then show how internal parameters contribute to 
the school’s cohesiveness. Furthermore, our results suggest a very interesting 
correlation between the degree of cohesiveness and the four avoidance patterns. 

The outline of this paper is as follows. Section [5] gives model description. 


We first recall the SDE model for fish schooling introduced in [26j, then newly 
inoculate a mechanism for obstacle avoiding into it. Section [3] presents four ob- 
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Stacie avoiding patterns. We thereafter investigate how these patterns change 
as the modeling parameters are tuned. Section H] explores fish school cohesive¬ 
ness. A scientific definition and measurement of cohesiveness are introduced. 
The relationship between avoidance patterns and school cohesiveness is then 
investigated. The paper concludes with some discussions of Section [5l 


2. Model description 


In 


26l| . we introduced a SDE model of the form 


dxi{t) 


= Vidt + aidwi{t), 



+Fi(x,,v,)'^dt, 



- Xj) 


(v* - Vj) 




( 1 ) 


for an A-fish system moving in the space {d = 2,3). Here, Xi{t) and Vi(<) 
denote position and velocity of the i-th fish at time t, respectively. And || ■ || 
denote the Euclidean norm of a vector, hence ||xi — Xj|| represents the distance 
between the i-th and the j-th fish. 

We regarded each fish as a moving particle in The first equation is a 
stochastic equation for the unknown Xj(t), where (Jidwi denotes noise resulting 
from the imperfectness of information-gathering and action of the z-th fish. In 
fact, = 1,2,..., N) are independent d-dimensional Brownian motions on 

some probability space. 

The second equation is a deterministic equation for Vi(t), where 1 < p < 
q < oo are fixed exponents; a and j3 are positive coefficients of attraction and 
velocity matching among fish, respectively. And r > 0 is a fixed number. If 
||xj — Xjll > r then the z-th fish moves toward the j-th. To the contrary, if 
||xi — Xjjl < r then the z-th fish acts in order to avoid collision with the j-th, r 
being thereby a critical distance (for details, see H)- 

Velocity matching of the z-th fish to the j-th also has a similar weight depend¬ 
ing on the distance ||xi — Xj jj. Degree of matching is higher when ||xi — Xjj| < r 
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for urgent reaction to avoid collision. When p is large, the attraction range 
is short. Meanwhile, when p is small, fish have ability to attract each other 
against a long distance. The exponent p thereby relates to a degree of how far 
the attraction reaches. Finally, the function Fi(xi,Vi) stands for an external 
force function acting on the Tth fish. 

The system © was studied analytically and numerically in the subsequent 
paper 171 of Quantitative arguments on fish schooling were also developed. 
Main advantages of using SDE models like ([T]) over others may be the conve¬ 
nience of mathematical technique. One can utilize the well developed theory 


of SDEs including numerical methods (see [l^). Its flexibility may be another 
advantage. As seen below, it is very easy to modify the model in the existence of 
an obstacle simply by introducing a suitable external force functions Fi(xi, v^) 
in (P). 

We now want to study from the viewpoint of obstacle avoiding. For this 
purpose, let us put a global obstacle with central point xq and radius p > 0 in 
and assume that fish move in the domain r2 = {xGK‘^:||x — xqH > p}. 
The surface of the obstacle is denoted by S = {x € : ||x — xc|| = p}. In 

addition to the rules (a)-(c) in Introduction, we introduce a behavioral rule in 
order to avoid collision with the obstacle: 


(d) Each fish executes an action for avoiding the obstacle according to a re¬ 
flection law of velocity with distance depending weight. 

On the basis of this rule, we are able to formulate the external force function 
Fi(xi, Vi) which influences the motion of the i-th individual. For a given (xi, Vi), 
where ||xi — xc|| > p, let li be a ray with origin x^ and direction v^, i.e., 

/i = {x e : X = Xi -I- svi, 0 < s < oo}. 

When li intersects S', we define Rf(xi, v^) as the reflection vector of at 
with respect to the tangential plane of S at y^, where y^ is the first intersection 
of li to S. In fact, is a vector whose opposite equals to the symmetric vector 
of Vi over the line connecting y^ and xq on the plane generated by and xq. 
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When li does not intersect S (including the special case when = 0), we put 
Rf(xi,v^) = Vi. Fig. 1 illustrates (x^,v^) and Rf(x^,Vi) in two-dimensional 
space. 

\ 

\ 

\ 


\ 

\ 

\ 



X. 

J 


Figure 1: Reflection vector Rf(xi,Vi) = 


Analogously to the velocity matching, we formulate 


Fj(xi, Vi) = -7 


RP RQ 

|xi -y»|l^ ||x, -y,||<3 


[vi - Rf(xj, Vi)], (2) 


where 1 < P < Q < oo are exponents, i? is a fixed distance, and 7 > 0 is a 
constant. 


If jjxi —yi II < R, then the z-th fish promptly reacts for matching its velocity to 
the reflection vector Rf(xi, Vi) to avoid collision with the obstacle. Meanwhile, 
if ||xi — yi|| > P, then the reaction to avoid the obstacle is less strong. If the 
ray li does not meet S, the individual takes no reaction with the obstacle. 
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The model for fish schooling in obstacle domain is therefore of the form © 
with the external forces Fi(xi, v^) (* = 1, 2,..., TV) given by ([2|). 


3. Obstacle avoiding patterns 

In this section, we observe four avoidance patterns of fish school based on 
our model o, where the external force functions v^) (T = 1,... ,N) are 

defined by ([5]). The effect of control parameters to these patterns is also inves¬ 
tigated. 


3.1. e-Graph and e, 9-Schooling 


Let us review notions of e-Graph and e, 0-Schoohng that were introduced in 
the paper [17 1. 


Definition 3.1 (e-Graph). Let Ki,\-i (i = denote a solution to 

SDEs (HI) . Let e > 0 be a fixed length. Regard the positions x^ (t) of fish at each 
time t as the vertices of a graph. Two vertices :}Ci(t) and Xj(t) are said to be 
connected by an edge of graph if and only if ||xi(t) — Xj(t)|| ^ e. Such a graph 
is called the e-graph of group at time t and is denoted by G£(t). 


Denote by n^ft) the number of connected components of Gg{t). When 
ng{t) = 1, the individuals form a single group. Meanwhile, when Ugit) ^ 2, 
they form ngft) sub-groups. 

Denote by crV{t) the variation of velocities: 


crV(t) 




1 

N 




2=1 


vWII^ 


0 < t < oo. 


N 

where v(t) = ^ is the average of all velocities of fish at time t. 

i=l 

Definition 3.2 (e, 0-Schooling). Let e > 0 and 9 > 0 be given. If a solution 
(x,,Vi) (* = 1, 2,..., A^) to dD) satisfies ns{t) = 1 and crY{t) ^ 9 for all suffi¬ 
ciently large t, say all t ^ T with some fixed time T > 0, then the state of fish 
is said to be in e, 9-schooling. 
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When the distance e > 0 and the tolerance 0 > 0 are specified, we simply 
say “in schooling” instead of “in e, 0 -schooling”. 

3.2. Obstacle avoiding patterns 

Set parameters as d = 2, = 20, a = /3 = 7 = 1. The exponent p is tuned 

from 2 to 4 keeping always the relation q = p + 1 and di = 0 for all i. The 
critical distance r is set by r = 0.5, and the radius of the obstacle is p = 1.2. In 
addition, e = 0.5 and 9 = 10“®. 

By performing preliminary computations, we first set a stationary state 
which is in e, 0-schooling. This state is set as the initial position of our compu¬ 
tations. The distance from the center of the school to the center of the obstacle 
is 3.5, and the line connecting these two centers coincides with the horizontal 
axis. The initial velocities are Vi(0) = (1.75, 0) for every i. Parameters for ob¬ 
stacle avoidance are set a.s P = p, Q = q, R = r. The school is thereby oriented 
toward the obstacle and will strike on it after a while. 

Fig. [5] shows the numerical results for p = 2, 3, 3.62, 4. Four different kinds 
of avoiding patterns are found. We will call them, Rebound, Pullback, Pass and 
Reunion, and Separation, respectively. 

Let us describe these four patterns of fish schooling. 

Pattern I (Rebound): The fish keep schooling throughout the obstacle 
avoiding process and the school rebounds off the obstacle. In order to keep being 
in schooling, they change their directions after the school touch the obstacle. 

Pattern II (Pullback): The individuals are once separated while ap¬ 
proaching the obstacle and stay around the surface of obstacle for a while. 
They then pull back off the obstacle to reform a school structure. 

Pattern III (Pass and Reunion): The fish pass the obstacle by spiffing 
to move along the obstacle surface. After passing it they reunite into a single 
school. 

Pattern IV (Separation): It is similar to Pattern III. But, after passing 
the obstacle, the subgroups have their own direction. 
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Figure 2: Obstacle avoiding patterns in two-dimensional space: Rebound, Pullback, Pass and 
Reunion, and Separation 


Similar results are observed in the three-dimensional case, i.e., d = 3 too. 
Numerical solutions to m show that there are similarly four avoidance patterns 
(I)-(IV). 

Fig. [3] represents four behavioral patterns of fish school while avoiding a 
static sphere obstacle in three-dimensional space. Fonr rows in this figure 
illustrate the time evolution of school starting from different initial speeds 
Vo = 0.3, 2.5, 6,13.5, respectively, while the other parameters are kept being 
constant in such a way that p = 2, a = /3 = 'y = 2, r=l. Initial positions of 
fish are set in s, 0-schooling, where e = 1 and 0 = 10“®. 
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Figure 3: Obstacle avoiding patterns in three-dimensional space: Rebound, Pullback, Pass 
and Reunion, and Separation 


3.3. Exponents p,q and obstacle avoiding patterns 

Set the parameters: d = 2, N = 20, a = (3 = 'y = 1, r = 0.5, p = 1.2, s = 0.5 
and 9 = 10“® as Subsection 13.21 The initial data (velocity and position) are 
also the same. 

The exponent p is however finely tuned from 1.001 to 8 with increment 
Ap = 0.001 and the relation g = p + 1 is always kept. In addition, P = p, 
Q = q, R = r. 

Our numerical results are presented in Table 1. This table highlights that 
there are critical values of p at which the type of avoidance patterns changes 
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from I to II, from II to III, and from III to IV, respectively. 


Table 1: Exponent p and patterns 


p [1.001, 2.100] [2.101, 3.371] 

[3.372,3.497] 

[3.498,8.000] 

Pattern I II 

III 

IV 


3.4- Speeds of school and obstacle avoiding patterns 

We set a fish group being in schooling. We investigate relations between the 
speed of school ||h|| and the type of performed pattern. 

Table [2] shows our numerical results in two-dimensional case. Here N = 20, 
Qf = /3 = 7=1, p = 2, r = 0.5, and the distance from the center of fish 
in schooling to the center of the obstacle is 3.5. The increment A||'(;o|| is also 
0 . 001 . 


Table 2: Speeds of school and patterns 


lluoll [0.001, 1.199] 

[1.200, 2.589] 

[2.590,4.866] 

[4.867,20] 

Pattern I 

II 

III 

IV 


Of course, if the initial speed is too large, then collision may happen because 
the large velocity makes fish have no enough time to adjust to avoid collision 
with other fish or with the obstacle. We do not consider this case. It is, however, 
interesting to know that solutions to o, where Fi{i = 1,2,..., N) are given by 
@, can blow up in a finite time. 

3.5. Critical distance r and obstacle avoiding patterns 

Let tune the critical distance r from 0.2 to 2.8 with increment Ar = 0.1. 
Other parameters are d = 2, N = 20, a = f) = ^ = 1, p = 1.2, P = p = 3, 
Q = q = 4,e = R = r, 6» = 10"®. 

In order to set the initial positions, we perform preliminary computations 
for each r in the free space where Fi = — 5vi. These computations provide 
stationary states in e, d-schooling for each r. 


12 








0, 


As pointed out by [17|, the geometrical diameter S 


N 


6 = max llxi — x|| 

l<i<N 


with X = ^i)/N 


of e, 0-schooling depends on r. It is thus natural to choose different radius of 
obstacle depending on r. In the present computations, the radius of obstacle 
is set as p = 2S{r)/3, where d{r) is the school diameter corresponding to r. 
The distance from the center of school to the center of obstacle is 8, and the 
line connecting these two centers coincides with the horizontal axis. The initial 
velocities are Vi(0) = (4,0) for all i. 

Our numerical results given in Table [3] show that as r increases, the type of 
patterns changes from IV to III, from III to II, and from II to I. 


Table 3: Critical distance and patterns 


r [0.2, 0.3] 

[0.4, 0.5] 

[0.6,2.0] 

[2.1, 2.8] 

Pattern IV 

III 

II 

I 


4. School cohesiveness 

In this section, we want to introduce a scientific definition of cohesiveness 
possessed by fish school as a nature of school. We then investigate the re¬ 
lationships between school cohesiveness and obstacle avoiding patterns. These 
relationships suggest that avoidance patterns can be used to visualize the school 
cohesiveness. 

4-1. Definition and measurement of school cohesiveness 

Let us first introduce a scientific definition of cohesiveness for the system o 
of SDEs in free space. 

Definition 4.1. School cohesiveness is the ability of group offish to form and 
maintain the e, 9-schooling structure against the noise imposed on the school. 
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In other words, how far the group can keep on e, 9-schooling as the magnitudes 
of the noises increase. 

This definition is given in a quantitative form. When e and 9 are specified, 
it is possible to quantitatively measure the cohesiveness of a group. 

Let us next give some examples of measuring cohesiveness of fish school by 
numerical methods. 

Consider a group of 50 fish moving in The parameters are set as p = 4, 
g = 5, Of = 4, /3 = 1. The external force functions are taken as Ti = —for 
i = 1, 2,..., iV. Initial positions Xi(0) are randomly located in a suitably small 
domain with null initial velocities Vi(0) = 0. The magnitude Ui = a \s a control 
parameter of simulation. 

We pick out 20 different trajectories of the Wiener process. For each value 
a, numerical computations for the solution Xi(t) and Vi(t) are performed in 20 
trials corresponding to these trajectories. 

Set e = r = 0.5 and 9 = 0.05. It is examined whether or not the states 
(xi(t),Vi(t))(i = 1,2,...,A^) are in £,0-schooling by fixing T = 30. In other 
words, it is checked that whether or not ne{t) = I and aV{t) < 9 for every 
t > 30. 

Starting with sufficiently small a, a is then increased with increment step 
0.001. When the e, 0-schooling structure is broken down at least for one sample 
trajectory of the Wiener process, the fish are considered to have lost ability of 
schooling. The critical value of a which is the largest value of a such that the 
group is still in e, 0-schooling can then be found. 

Fig. 0] shows that when a = 0.02, the group builds up the e, 0-schooling 
structure. Here, the sub-figure on the left shows the e-graph of the group at 
time t = 30, the positions Xi(t) being drawn by dots, and the edges of graph 
by lines. The right hand side sub-figure illustrates the variance of velocity as a 
function of t, where the horizontal line represents the level 0 = 0.05. 

The group loses schooling ability when a = 0.06. In fact. Fig. [5] demonstrates 
that ne{t) = 2 for t = 30. Schooling ability of the group is also lost when 
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a — 0.069. Fig. ini highlights that crVS(t) > 0 for t = 30. 

By these methods, we can finally compute the critical value a = 0.051. The 
group is in e, 0-schooling whenever a < a, whereas it loses schooling ability 
whenever a > a. 



Figure 4: e-graph and o-VS(t) for a = 0.02 


4-2. Relations between exponents p,q and school cohesiveness 

Let tune p as p = A, 3.62, 3, 2, keeping the relation q = p + A and other 
parameters as before. Our computations show that the numerical critical value 
of (T, namely the cohesiveness, changes as ir = 0.051, 0.055, 0.056, 0.063, respec¬ 
tively. 

The exponent p, as explained in Section [21 shows a degree of range how far 
the attraction is effective. It is therefore very natural that as p decreases, the 
attraction range extends and enhances the cohesiveness of group. 

4-3. Relations between critical distance and school cohesiveness 

Let us tune r from 0.5, 0.6 and 0.7 by taking e = r, p = 4, g = 5, 0 = 0.05, 
T = 30. The corresponding cohesiveness is found as tr = 0.051,0.053,0.054, 
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Figure 5: e-graph for cr = 0.06 



Figure 6: (TVS(f) for cr = 0.069 

respectively. This means that cohesiveness is enhanced as the critical distance 
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increases. 


4.4- Relationship between school cohesiveness and obstacle avoiding patterns 

In Subsections 13.21 and 13.31 we already show that suitable tuning of the 
exponents p and q = p + 1 provides different obstacle avoiding patterns. It is 
possible to interpret this fact as follows. 

Note that p is concerned with the range of attraction among fish. If p is small, 
then the attractive force reaches a wide range beyond the critical distance r. In 
contrast, if p is large, then the attraction is only available in a neighborhood of 
the disk of radius r. In Subsection we verify that when other parameters 
are fixed, the cohesiveness of school increases as p decreases. 

In the numerical examples in Subsection 13.21 when p = 2, the school has 
very strong cohesiveness and rebounds off the obstacle. When p = 3, the school 
still has strong cohesiveness and can keep schooling. But the fish are spread 
on the surface. When p = 3.62, the school cohesiveness becomes smaller. The 
fish can no longer keep being in schooling but it is strong enough to reunite the 
members into a school. When p = 4, the school cannot keep being in schooling 
and is separated into two clusters after passing the obstacle. 

If these interpretations are reasonable, then the four obstacle avoiding pat¬ 
terns can be used to measure the cohesiveness of school approximately. For 
example, the cohesiveness can be easily categorized into four classes. 

5. Conclusions 

Obstacle avoiding patterns and cohesiveness of fish school have been studied. 
In fact, our mathematical model for fish schooling in the space with obstacles 
provided four clear avoidance patterns: Rebound, Pullback, Pass and Reunion, 
and Separation. The shift from one pattern to another due to the change of 
parameters has also been investigated. 

Furthermore, our definition for school cohesiveness suggested that the cohe¬ 
siveness can roughly be measured by using the four types of patterns. In other 
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words, we could visualize the “strength” of the school by connection with the 
patterns. Quantitative understanding for the cohesiveness of swarming systems 
may in turn provide useful information in designing artificial self-organizing 
systems or intelligent systems [s, 22|. 


Our attempt may however not be complete, since the emerging patterns 
depend not only on the internal parameters but also on the environmental con¬ 
ditions, for instance the radius of obstacle. If the geometrical diameter of school 
is relatively larger than the obstacle’s radius, then the school striking the ob¬ 
stacle may pass over and reunify. To the contrary, if the school’s diameter is 
relatively smaller than the radius, then the school may rebound off. 
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